Scaling behavior of the adiabatic Dicke Model 
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We analyze the quantum phase transition for a set of A''-two level systems interacting with a 
bosonic mode in the adiabatic regime. Through the Born-Oppenheimer approximation, we obtain 
the finite-size scaling expansion for many physical observables and, in particular, for the entangle- 
ment content of the system. 



^ . 

o ■ 
o ■ 



> 

cn 
o 

o 

■ 

a • 

^ : 



PACS numbers: 64.60.Fr, 03.65.Ud, 05.70.Jk, 73.43.Nq 

Introduction. Many-body systems have become of 
central interest in the realm of quantum information the- 
ory as training grounds to study static, dynamical and 
sharing properties of quantum correlations. It has been 
natural, then, to employ the entanglement as a tool to 
analyze quantum phase transitions, one of the most strik- 
ing consequences of quantum correlations in many body 
systems [lj,y,|3|. In this letter, we study the Dicke super- 
radiant phase transition and obtain the finite size scaling 
behavior of the quantum correlations at criticality. 

The Dicke model (DM) describes the interaction of N 
two-level systems (qubits) with a single bosonic mode 
and has become a paradigmatic example of collec- 
tive quantum behavior. It exhibits a second-order phase 
transition 0, which has been studied extensively 
The continued interest in the DM arises from its broad 
application range ^ and from its rich dvnamics, dis- 
playing many non-classical features [ifl |lj [13, ll3 ■ The 
ground state entanglement of the DM has been recently 
analyzed flit . ,15^ .16] , and some aspect of its finite size 
behavior has been obtained numerically. Finally, Vidal 
and Dusuel, obtained the critical exponents by a 

modified Holstein-Primakoff approach. 

The exact treatment of the finite-size corrections to 
the Dicke transition is quite complicated and the study 
of some limiting cases can be useful. In this Letter we 
analyze the case of N qubits coupled to a slow oscilla- 
tor. We discuss on equal footings both the finite size and 
the thermodynamic limit of the model, thus allowing to 
obtain the phase transition as well as its precursors at fi- 
nite N. Indeed, we obtain the dominant scaling behavior 
for some entanglement measures and the entire 1/N ex- 
pansion for all of the relevant physical observables, such 
as the order parameter. Concerning the quantum cor- 
relations, we argue that what is really relevant for the 
critical behavior is the bi-partite entanglement between 
the oscillator mode and the set of two level systems. 

Adiabatic limit. The interaction of TV identical qubits 
with a bosonic mode is described by the Hamiltonian 

(1) 



+ AS^ + A(at + a)S,/VN 

where uj is the frequency of the oscillator, A is the transi- 
tion frequency of the qubit and A is the coupling strength. 



The 5"s are the total spin observables, S^, = J2iLi'^l' 
where cr* is the v-th the Pauli matrix for the z-th qubit. 
H is equivalent to the Dicke Hamiltonian 4], obtained 
after the rotation e"'^"/*. 

We assume a slow oscillator and work in the regime 
A ^ w by employing the Born-Oppcnhcimer approxima- 
tion. As detailed in 18], the procedure can be followed 
more plainly by rewriting the Hamiltonian of Eq. |^ as 



^^2 



DS., 



LQ 



N 



(2) 



where we have introduced the dimensionless parameters 
D ~ 2A/(jj, and L — 2^/2X/uj, together with the oscilla- 
tor coordinates Q = {a'^ + a)/ V2, and P = i{a^ — a) / \/2. 

The basic assumption of the well-known adiabatic ap- 
proximation is that the state of a composite system with 
one fast and one slowly changing part can be written as: 



(3) 



\x{Q)) is the eigenstate of the "adiabatic" qubit equation 
for each fixed value of the slow variable Q, 



Ds. + ^s,]\xm = Emxm, (4) 



and can be written as the direct product of the eigen- 
states of each single qubit 

\xiQ)) = lx(Q))i ® \x{Q))2 ® • • • ® \xmN . (5) 

As the qubits are identical, the lowest eigenstate of Eq. 
has the form 
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ixo(o)) = { ^ [A-m+) A+m-)] } , (6) 



where |±) are the ±1 eigenstates of a^, and 

A±(g) = ./i± 



LQ 



NQiQ) 



(7) 
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The eigenvalue corresponding to this state is given by 



EoiQ) = -NQiQ) 



L2Q2 
N 



(8) 



This energy eigenvalue contributes an effective adiabatic 
potential felt by the slow bosonic mode. The total po- 
tential for the variable Q is, therefore 



UoiQ) ^ [Q' ~ NSiQ)] 



(9) 



Introducing the dimensionless parameter a — L^/2D, 
one can show that for a < 1, the potential Uo{Q) can 
be viewed as a broadened harmonic potential well with 
its minimum at Q = 0. For a > 1, on the other hand, 
the coupling with the qubit produces a symmetric double 
well with minima at Q = ±Qo = ±^/N D\/ a'^ — 1/L. 

Ground state properties. In order to obtain the fun- 
damental level of the coupled system, the last step in the 
adiabatic procedure is the evaluation of the ground state 
wave function for the oscillator, (j)o{Q), to be inserted in 
Eq. . This wave function satisfies the one-dimensional 
time independent Schrodinger equation 



uj d 
'2dQ 



MQ) = £oMQ), (10) 



where eq is the lowest eigenvalue. 

The ground-state properties can be easily studied in 
the adiabatic limit. It turns out that all of the expec- 
tation values can be expressed in terms of the quantity: 



00 (Q) 1 



2a 



dQ, 



(11) 



which is shown below to depend only on a and ND. 

For the average values of the various components of 
the total spin, one gets 



N 



(5,) = (5.)=0, 
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(SI) (S'y) 



1 
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(12) 



(13) 



(14) 



Furthermore, the order parameter can be obtained from 
(H) Lu/(P^ + Q^ 



N 



N 



- L>$i(iV) 



(15) 



In the thermodynamic limit (A^ — > oo) one gets the 
well-known second-order quantum phase transition at 
a = Q!r = 1, for which: 



N 



-1 (a<l) 

-i («>i), 



(16) 




FIG. 1: The behavior of {Sx)/N and of the ground state 
energy (inset) as a function of the parameter a, for D = 10. 
Different curves correspond to different values of A'^. 
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and, finally 



-D {a < 1) 

(o>i), 



-0, 
(17) 



(18) 



(19) 



where Eq ~ 2eo /uj. Numerical results for the ground 
state energy and the x-magnetization are plotted in 
Fig.© as a function of the parameter a, for different 
values of N, in comparison with the results for N ^ oo. 

In what follows, we will obtain an analytic expression 
for the N dependence of ^i,{ND), and, consequently, the 
scaling relations for all the observable introduced above. 

Finite-size scaling exponents at the critical point. In 
the adiabatic regime, i.e. for large D, the Schrodinger 
equation Hl()|l can be approximately rewritten as 



2ND 

eoia,ND)MQ;a,ND) (20) 



with eo(a, ND) = Eo{a, ND) + ND, ll9|. 

Even thought the ground state energy and all of the 
physical observables appear to depend on the two pa- 
rameters a and ND, this problem can be simplified in a 
single-parameter one with the help of Symanzik scaling 
|20| , after transforming Eq. H2UI) into the equivalent form 



dq^ 



</'o(<z;C) = eo(C)0o(g;C) (21) 



where q = Q (2^)'^', while C = («c - «) is 

the only remaining scale parameter. 



3 



For the ground state energy one obtains the relation 

Eoia,ND)^-ND+i^-^j eo (C) • (22) 



Taking the hmit a ^ Oc (or ^ — > 0), Eq. (|21|l becomes 

+ q*^Mq\0)^eo{0)Mq;0) (23) 



whose lowest eigenvalue is found to be eo(0) ~ 1.06036. 
For C 7^ (but near the critical point) we can resort to 
perturbation theory and obtain the ground state energy 
as an expansion in powers of C, 



eo 



(c) = E/^»c". 



(24) 



It is easy to show that /3o — eo(0) and /3i = 
^^^q^(t>l{q;Q)dq = €'^{0) ~ 0.36203. We demonstrate 
below that these /3's enter not only the energy, but also 
the finite N expansion of every physical observable. 

In order to compute the average values listed in Eqs. 
(|12I15|) at the critical point, we expand the integral 
^v{ND), given in Eq.|ni), as 



M7Vi^)^i + ^(g^) 



2a'^{v- l)v 



{Q') (25) 



Taking the system energy -Eo(a, ND) from Eq. H22|) . we 
exploit the Feynman-Hellman theorem to obtain 

^^^%^^-(0^) + ^(Q^), (26) 



dEQ{a,ND) 



-1 



d{ND) ^ 2N^D^ 

Using Eq. (|24|l . we get the critical point values 

00 



(27) 



(Q2) = {2ND)^/^ , (g4) = ^ {2NDf'^ (28) 



^,{ND) ~ 1 
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{2NDfl^'^ 3{2ND)^/ 
Substituting into Eqs. l(T^ - l(0|l . we have 
(S.) , , 2/3i 2/3o 



N 

JV2 



{2ND)y^ {2NDY/^ 



N -1 
N 



4/3i 



16 



/3o 



(2iVL>)2/3 3 (27Vi:>)4/3 



Finally, the order parameter is given by 



ND 



2/3i 



/3o 



(2iVi:')2/3 3(2iVi:')4/3 



(29) 

(30) 
(31) 

(32) 



In Fig. |(2Jl we make a comparison with the results ob- 
tained from the numerical solution of the Schrodinger 
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FIG. 2: Scaling of the ground state energy and spin average 
as a function of A'^ at the critical point a = 1, for D — 10. 



equation (|10|l . One can see that the leading finite size- 
corrections for eo/ND and {Sx)/N scale indeed as N~^/^ 
and N~'^^^, respectively. When the system size is too 
small, the next to leading orders in these expansions be- 
come important and this explains the discrepancy with 
Ref. 01 , already pointed out in 0| . 
The knowledge of /3o and of the first two moments is 
sufhcient to determine recursively all the others {k > 2): 



{Q 



k+4\ 



{2NDY/^ k + i 



4(fc + 3) 



■{2NDYi-\Q 



l/3/r)A;+2\ 



(33) 

This relation allows one to compute higher order terms 
of the finite size expansion. 

We conclude this part by pointing out that the os- 
cillator state (}io could also be found by a different re- 
scaling procedure, i.e. by rewriting the Hamiltonian as 
H{^) = p"^ +q^ + ^i^jj ^ ^ a^/2ND{a^ - af/^. 
In this case, however, the perturbation expansion would 
yield a power series that diverges very strongly for every 
^ 7^ 0. This would be the equivalent in our approach of 
the method of Ref. 0| , where the scaling exponents are 
obtained by arguing that there can be no singularity in 
any physical quantity at finite-size. Our leading order 
results agree with those reported in Ref.|3 except for 
the exponent of {Sy)/N'^, on which we comment below. 

Entanglement. To start the discussion on quantum 
correlations in this model, we point out the peculiar na- 
ture of the adiabatic state of Eq. namely, the fact 
that, once the oscillator is traced out, the TV qubits re- 
main in a purely statistical mixture, without the presence 
of any entanglement, neither of pairwise, nor of multi- 
partite nature. In particular, the concurrence between 
any two qubits is zero. This is due to the adiabatic hy- 
pothesis 3> 1, which strongly suppresses the energy 
exchange between qubits, mediated by the oscillator. 

Since we obtain all the relevant features of the Dicke 
transition despite the fact that concurrence is neglected, 
we can conclude that its presence is not really essential 
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FIG. 3: The tangle rjv between the oscillator and the A'^ 
qubits as a function of a, for D = 10 and different values 
of A'^. The inset shows the scaling of rjv with TV at a = 
Oc = 1, compared with the analytic expression of Eq. I38II 
(continuous plot). 



to describe the large N behavior. In fact, as reported 
by Ref. 0|, some degree of pairwise entanglement is 
present even in the regime D ^ 1 and therefore one 
should expect that the adiabatic approach fails in de- 
scribing pairwise correlations. Interestingly enough, this 
is not entirely the case as, for large N, the only two-point 
correlation function for which we obtain a different scal- 
ing behavior is (Sy), and this difference is just enough to 
set the concurrence to zero. 

As argued above, the quantum correlations that really 
matter for the phase transition are those involving the os- 
cillator. Below, we evaluate i) the entanglement of each 
qubit with the rest of the system (as there is no entan- 
glement among qubits, this is due to correlations with 
the boson mode), and ii) the amount of entanglement 
between the oscillator and the entire set of qubits. 
After tracing out all of the other degrees of freedom, each 
qubit is found in the same state pi, and participates of 
the entanglement (as measured by the tangle) 



Tl 



1 - TtpI 



1 



iV2 



(34) 



In the thermodynamic limit, this is zero in the normal 
phase, while one gets ti = 1 — ^ for a > 1. For large 
but finite N, ti is non-zero even for a < 1; but at the 
critical point it scales as ti ~ (^2ND)^/'-^ • 

We adopt the linear entropy also to quantify the entan- 
glement between the oscillator and the N qubits. Thus 



r/(l -Tr{p^ 



(35) 



where pN — Trosc|V')(V'l is the reduced qubit state, while 

to bound to 



the pre- fac t or is chosen to he rj — ^rr: 
1, see Ref. |^. Explicitly, we have 



Trp; 



dQdQ' ^1{Q)4{Q') 



e(g)e(Q') 



N 



(36) 



which can be expressed as a power series in 

Fig. (PJ shows Tat both for finite N and for N 
In the thermodynamic limit, we find 



1 



1 - 



1 



(a<l) 
(a>l), 



(37) 



1. 



which shows a cusp at the critical point, where Tq, 
For finite N , this singular behavior is rounded and r is 
quenched. When N is large, the entanglement scales as 



Tat (a = ttc) ~ 1 - 



(38) 



(2D)i/3iVi/6 ' 

where K — J dq(t>^{q\ 0) ~ 0.46, and 4)i){q\ 0) is the nor- 
malized solution of Eq. (|23|) . This result is shown in Fig. 
(O to agree with the numerical evaluation of Eq. H35() . 

To summarize, we have described the finite size scal- 
ing of quantum correlations in the Dicke model, obtaining 
the finite N behavior of the spin components, of the order 
parameter and of the ground state entanglement in the 
adiabatic regime D ^ 1. We also discussed the crucial 
role of the entanglement involving the oscillator at crit- 
icality, giving its expression both in the thermodynamic 
limit and for a finite size system. 
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